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ABSTRACT 

We present spatially resolved kinematics of the interstellar Nai D AA5890,5896 doublet 
absorption and 12 CO (1— >0) emission across the inner ^2 x 1 kpc of the disk of M82. 
These data were obtained with the DensePak IFU on the WIYN telescope and the 
Caltech Owens Valley Radio Observatory (OVRO) millimetre array. By measuring 
the Nai and CO (and Ha kinematics from a previous study) at the same spatial 
resolution, and employing the same line fitting method, we have been able to make 
meaningful comparisons between the ionized, neutral and molecular gas phases. We 
detect a component of the Na I line throughout the inner disk with velocities that are 
forbidden by the known galactic rotation. We interpret this as originating in counter- 
rotating or perhaps inflowing material. In the southern plume, we find clear evidence 
of entrained CO gas with kinematics consistent with that of Ha. On the northern 
side, the CO kinematics appear to trace more static clouds in the inner halo that 
could be pre-existing or tidal in origin. We find no evidence that Nai absorption is 
kinematically associated with the outflow. We conclude that a combination of lack 
of velocity resolution and confusion of due to the high inclination of the system is 
acting to prevent detection. Thus, in the search for neutral outflows from galaxies, 
the signature high velocity components may easily be missed in observations at low 
spectral resolution and/or sensitivity, and particularly so in highly inclined systems. 

Key words: galaxies: individual (M82) - galaxies: starburst - galaxies: ISM - ISM: 
kinematics and dynamics. 



1 INTRODUCTION 

M82 is the archetype nearby (3.6 Mpc, 1" 



17.5 pc; Freed 



man et al. 19941 starburst galaxy (O'Connell & Mangano 
1978||Q'Connell et al.||1995| ). The current (~10 Myr) star- 
burst activity is concentrated in a ~500 pc (~30") region 
centred on the nucleus. From Hubble Space Telescope (HST) 
imaging, the starburst is known to consist of many young 
massive star clusters ( O'Connell et al.(T9 95 ; McCrady et al. 
2003 1 , and it is the energy from these clusters that drives the 
famous Ha- and X-ray-bright, kpc-scale, superwind (Shop- 



bell & Bland-Hawthorn 



ct al. 2003; Strickland & 
2007b 



19981 lOhyama et al.||2002| |Stevens 



|2009a). 
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The bipolar outflow ( Shopbell & Bland-Hawthorn 1998 



Engelbracht et al. 2006 



Mutchler et al. 2007) can be ob- 



served easily due to the galaxy's almost edge-on inclination 
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(i ~ 80°; |Lynds fc Sandagej [1963] |McKeith et al.||1995| ). 
McKeith et al. (1995 hereafter |McK95 1 presented deep op- 



tical and near-infrared long-slit observations of emission and 
absorption lines along M82's minor axis, and found clear 
line splitting in their emission lines. This was interpreted as 
originating from the front and back walls of a cone-shaped 
structure, inclined such that the southern cone is directed 
toward the observer. 

Radial velocities along the major axis have been stud- 



ied many times in the optical/near-IR (e.g. 


Burbidge et al. 


1964| |McKeith et al.| ] 


1993 Sofue et al. |1998 Shopbell & 


Bland-Hawthorn| 1998 


Castles et al. 1991 


Wcstmoquette 


et al.||2007b[), IR (e.g. 


Larkin et al. 1 1994 Achtermann & 


Lacy||1995|) and radio (e.g. Wills et al. 2000 Walter et al. 


2002| |Rodriguez-Rico et al. 20041. Overall, these studies 



dient, and an outer region >±25" of constant (or slowly 
increasing) velocities. These gradients are centred on the 
systemic velocity of the galaxy, v sys — +200 kms~ . 
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The steep inner region is attributed to the a;2-orbits of 
a ~1 kpc long stellar bar, known from near-infrared and Hi 
studies I Larkin et al.|1994 Achtermann fc Lacy|1995 Wills 



et al.|2000 |. The interaction between the outer x\ and inner 
X2 orbit families has resulted in the build up of a torus of 
molecular and ionized gas and dust at a radius of ~250 pc 



(Achtermann & Lacy 1995 


Weifi et al.|2001 1. 




Westmoquette et al. 


(2009a hereafter 


Paper 1 1 



ined the ionized gaseous and stellar dynamics of M82's disk 
and inner wind via DensePak and Gemini/GMOS spatially 
resolved optical integral field unit (IFU) spectroscopy. We 
found clear evidence for the wind emanating from the entire 
starburst zone. We also found a broad component in the ion- 
ized emission lines that we interpreted as tracing turbulent 
mixing interfaces between the hot wind and photoionized gas 
phases. In Westmoquette et al. (2009b, hereafter Paper if] ) 



we examined the nebular properties in more detail, and de- 
scribed the discovery of a distinct channel > 100 pc in length 
identified within the inner wind, traced by the ionized gas 
and characterised by its dynamics and density morphology. 
The central starburst is believed to be fuelled by the 



(Aft, 



2.3 



x 10 s M ; 


Weifi et al. 


2001 


Garcfa-Burillo 


Walter et al. (|2002), however, discovered sig- 



nificant amounts of molecular material (1.3 x 10 Mq) out- 
side the central 1 kpc disk. Some streamers were found to 
be well correlated with optical absorption and Hi features, 
implying that the molecular gas has been partly disrupted 
by the tidal interaction with M81. However, CO line split- 
ting was also detected for the first time, implying that there 
is also significant molecular material entrained in the wind 
(~3 x 10 s Mq). 

While much work has been done on the distribution 



of ionized gas in nearby starburst outflows (e.g. 


Marlowe 


et al.|1995; Martin 


1998 


|Westmoquette et al.|2007a||Sharp 


& Bland-Hawthorn 


2010 


| Westmoquette et al. 2011 1, in com- 



parison very little has been done (particularly on spatially 
extended scales) on the neutral phases in these same nearby 
examples. Cool neutral gas is most easily probed through 
absorption lines. Due to their advantage in (1) being able 
to be seen against a bright continuum source (continuum 
surface brightnesses are typically much larger than that of 
the emission lines), and (2) mitigating the ambiguity in the 
radial flow direction, absorption lines have been used ex- 
tensively to probe outflows in more distant sources (z > 1, 



e.g.|Shapley et al.| 


2003| |Fox et al.|2007||Weiner et al. 2009 


Rubin et al.|2010 


Kornei et al.|2012 


Martin et al.|2012 


1. 



Since, in these samples, the absorption kinematics are 
often treated as direct probes of the outflow speed, it is im- 
portant to examine whether this really is the case in nearby 



systems. Schwartz & Martin (20041 presented Keck/HIRES 
Nal D absorption line observations of a sample of nearby 
dwarf starbursts and found that the kinematics of Nal D 



and Ha were generally very different. Rupke et al. (20051 
compared the ionized emission and Nal D absorption line 
profiles for a sample of 17 starburst-dominated ULIRGs, 
and also found a complicated situation. Only in some cases 
blueshifted wings were detected in both phases. In half of 
those, the ionized gas had a higher maximum velocity, while 
in the other half they were equivalent. The situation, there- 
fore, is far from clear. 

In this paper we aim to determine the kinematics of 



both the cool neutral and molecular media across the en- 
tire inner ~2 x 1 kpc of M82's disk using spatially re- 
solved (3" resolution) observations of the interstellar Nal 
D AA5890,5896 doublet absorption and the 12 CO (l->0) 
(2.6 mm, 115 GHz) emission line. By comparing these re- 
sults to those of the ionized gas (made at the same resolu- 
tion; Westmoquette et al.|2009a |, we can make a direct and 
meaningful comparison between the three gas phases. 

2 OBSERVATIONS 

2.1 DensePak observations 

On 14th April 2001, we observed four fields covering the 
stellar disk of M82 with the DensePak instrument ( |Barden| 
et al.|1 998 ) , a small fibre-fed integral field array attached to 
the Nasmyth focus of the WIYN (Wisconsin, Indiana, Yale 
and NOAO) 3.5-m telescope. It has 91 circular fibres, each 
with a diameter of 300 um or 3" on the sky; the fibre-to- 
fibre spacing is 400 um making the overall dimensions of the 
array 30 x 45 arcsec. The format of the array, including the 
gaps between the fibres, is represented in the maps presented 
below. DensePak's fibre bundle is reformatted into a pseudo- 
slit to feed the WIYN bench-mounted echelle spectrograph. 

The spectrograph was set up to give a dispersion of 
0.46 Apix -1 over the wavelength range of 5820-6755A. This 
allowed us access to a number of optical nebular lines (anal- 



ysed in Paper I and Paper III and the Nal D AA5890,5896 



doublet which we examine here. Nal D absorption origi- 
nates from warm neutral gas in the interstellar medium 
and stellar atmospheres. We refer the reader to |Paper lj 
for further details regarding the IFU position coordinates, 
exposure times and data reduction methodology. Fig. [I] 
shows the position of the DensePak fields overlaid on an 
HST /ACS F814W (I-band ) image. The instrumental reso- 
lution is 44.7 ± 3.1 kms" 1 (IPaper 1 1. 



2.1.1 Decomposing the Nal line profiles 

Nal absorption is detected in every fibre of the DensePak 
fields, and the high S/N and spectral resolution of our data 
have allowed us to quantify the line profile shapes to a high 
degree of accuracy. Because of the coincidence of telluric 
Na I emission (mostly from the street lights in Tucson) with 
the interstellar absorption lines, accurately measuring the 
absorption profile shape is highly dependent on the sky sub- 
traction quality. Thus, to avoid the introduction of artefacts 
that are always associated with attempting to remove the 
sky lines, we here perform our analysis on the original spec- 
tra (i.e. with no sky subtraction). 

For each spectrum we simultaneously modelled the Na I 
doublet absorption together with the He I A5876 nebular 
emission line and the three telluric emission lines at 5888.05, 
5889.95, and 5895.92 A (the latter two being the Nal dou- 
blet) using using the IDL \ 2 fitting package PAN (Peak ANal- 
ysis |Dimeo|[2"005"l |Westmoquette et al.1|2007a|). Milky Way 



absorption was filled in by the telluric emission so was not 
included in the fit. In order to limit the number of free pa- 
rameters as much as possible, we fitted the Nal doublet 
(both emission and absorption) as a single model consist- 
ing of two Gaussians with wavelength separations equal to 
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the lab-measured difference, and equal line widths specified 
to be greater than the instrumental resolution. The fluxes 
of the two lines of the doublet were left unconstrained. We 
fitted each of the Na I absorption profiles with a single and 
double component initial guess. Multi-component fits were 
run several times with different initial guess configurations 
(widths and wavelengths) in order to account for the varied 
profile shapes, and the one with the lowest x 2 fit statistic 
was kept. However, we note that the x 2 minimisation routine 
employed by PAN is very robust with respect to the initial 
guess parameters. 

To determine whether the Nal doublet needed one or 
two components to fit its shape, we used a combination of 
visual inspection, the x 2 , and the fit residuals. In the major- 
ity of spectra, two components were needed. In these cases, 
we named the stronger absorption component CI and the 
weaker C2. Fig. [2] shows some example spectra in the Nal 
region, showing simultaneous fits to the He I A5876 emission 
line, the Na I absorption lines and the three telluric emission 
lines (note that He I is not detected in two of these exam- 
ples). These spectra were extracted from the spaxels labelled 
with the corresponding letters in Fig. [4] (see below). 

Due to the constraints imposed on the telluric line fits, 
formal errors on the Na I fit results are not straightforward to 
calculate. We estimate that uncertainties in the radial veloc- 
ities of the Nal results are 5-10 kms -1 for high S/N spaxels 
and remain ^30 kms -1 for the majority of spectra. FWHM 
errors range between ^10 kms -1 for high S/N spectra and 
remain below ~50 kms -1 for most spaxels. Uncertainties in 
the absorption equivalent widths (W cq ) of the components 
range between 0.1 A for the strongest absorption to ~0.5 A 
for the weaker components. In some cases, these high un- 
certainties arise from the A5896 (redder) Nal emission line 
falling on the red wing of the bluer (A5890) absorption line. 

Fig. [3] shows the Nal W oq maps for components CI 
and C2. The W cq of Nal C2 is a factor of ~3-5 times 
smaller than CI. As a check on the physicality of our fi- 
nal results, we measured the ratio of the equivalent widths 
of the two Nal lines for each fitted component, since a ra- 
tio of Weq,589o /W eq ,5896 = 1 represents the optically thick 
limit, whereas a ratio of 2 represents the optically thin limit. 
As expected, in almost every case the ratio is ~1. However, 
in spectra where the uncertainties in the W cq of the weaker 
component are high, this ratio is often <1 (Fig. [2]l is an 
obvious case). While this reflects the inherent difficulty of 
fitting such a complex and blended set of lines, we believe 
that, although the equivalent widths of the weaker compo- 
nents may be uncertain, the kinematics are accurate. 



2.2 OVRO CO map 

Observations of M82 in 12 CO (2.6mm, 115 GHz) were 

made with the Caltech Owens Valley Radio Observatory 
(OVRO) millimetre array over a a 2.' 5 x 3.'5 field at a spatial 

(120021. 



resolution of 3. 6, and were presented in Walter et al 



More details on the setup and data reduction methodology 
are given in that paper. Here we use these data to compare 
to our ionized and neutral phase measurements. 

First we cropped the cube to the region coincident with 
our DensePak data, then binned it spatially to a resolution 
of 3" to match that of the DensePak fibres. We then masked 
any spaxels with integrated fluxes below 15 Jy beam -1 



(equating to a line S/N below ~3) in order to ensure that ev- 
ery spatial element in our final cube contained a significant 
CO line detection. 

In order to produce measurements for the CO data 
that are comparable to our ionized and neutral gas mea- 
surements, we fitted the CO line using PAN in much the 
same way as we have done for the ionized emission lines JPa-l 



per 1 1 and Na I absorption lines. We found we required two 



Gaussian emission components in the majority of spaxels in 
the cropped and masked cube. Whenever two components 
were needed, we assigned the component whose radial veloc- 
ity was closest to v ByE (+200 kms -1 ) to CI and the higher 
velocity component to C2. We found that this resulted in 
a more coherent velocity map than if we had followed the 
convention adopted for Na I (where we assigned the stronger 
component to CI and the weaker to C2 regardless of veloc- 
ity). f\ Fig. [H] shows the resulting line centroid velocity maps 
for the two components, and Fig. [5] shows some example fit- 
ted CO line profiles extracted from selected regions labelled 
with the corresponding letters in Fig. [6] 



3 RESULTS 

3.1 Radial velocity maps 

In Fig. [4] we present the Na I line centroid radial velocity 
maps for the two identified line components. Recall that 
for Nal, the fainter component of any two-component fits 
was assigned to C2. The signature of galactic rotation is 
seen in the CI component, with approaching velocities in 
the west and receding velocities in the east, although the 
velocity gradient is much shallower than that of the stellar 
or ionized emission lines ( |Paper i] and see below) . In |Pa-| 
|per I| we confirmed a 12° offset to the position angle of the 
ionized gas rotation axis, compared to both the stellar ro- 
tation axis and photometric major axis, that had first been 
seen by Rodriguez- Rico et al. (20041. The line connecting 



the maxima and minima of this offset rotation axis is shown 
with a dashed line in the Nal CI map (Fig. H|. Despite the 
shallower velocity gradient compared to that of the ionized 
gas (see below), it is clear that the rotation seen in Nal 
CI is consistent with same offset axis (blueshifted velocities 
in the west peak below the major axis, whereas redshifted 
velocities in the east peak above the major axis). 

The dynamics of the fainter Nal component (C2) are 
much more complex. To the west of the minor axis, C2 is 
redshifted compared to CI (e.g. Fig. [2Jd). To the east, some 
spaxels contain a blueshifted C2 component whereas others 
show a redshifted one. In position 3, a handful of spaxels 
along the midplane exhibit a blueshifted C2 component (e.g. 
Fig. [5J;), although most C2 velocities in the south of the 
position are redshifted with respect to CI (e.g. Fig. [2ji). 

In Fig. [ts] we present the 12 CO emission line radial ve- 
locity maps. Two line components are detected throughout 
the central 80"x60" (1400x1050 pc). Unlike the flux-based 



1 In Nal, the stronger component dominates the absorption pro- 
file and is the one that mainly follows galactic rotation (see 
Sect |3.l) , whereas in CO the relative strength of the two com- 
ponents varies much more, meaning that a velocity-based assign- 
ment is more appropriate. 
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component assignments we used for Nal, for CO we assigned 
the components based on velocity, with the one closest to 
as CI and the higher velocity component as C2. The ro- 



( 1993 



tation of the central CO torus ( Shen fc Lo|1995 I can clearly 
be seen, with the redshifted receding side in CI and the 
approaching blueshifted side in C2. 

The differences between the Nal and 12 CO maps and 
the Ha maps from |Paper I| clearly indicates that the gas kine- 
matics within this central starburst region are very complex. 
To investigate these results in more detail we define a num- 
ber of pseudoslits along and parallel to the major and minor 
axes from which we extract the line component velocities. 
These are discussed in the following sections. 



3.2 Major axis position-velocity graph 

Position-velocity plots extracted from two 7"-wide pseu- 
doslits positioned along the galaxy major axis and offset 
15" to the south are shown in Fig. [7] (left and right panels 
respectively). The location of these pseudoslits are shown in 
both Fig.[4]and[6] We plot our new Nal and 12 CO measure- 
ments together with the Ha measurements from |Paper 1} 
the [Nn]A6583 velocities from Castles et al. (19911, and the 
[S in], P10, and Can velocities from McKeith et al 
+ symbols). 

A close inspection of these two plots reveals a number 
of points of interest. Between —20 to +5" on the major 
axis (the eastern side of the inner X2 bar region) the Ha, 
Nal and CO measurements are mostly in good agreement, 
indicating that we are seeing molecular, neutral and ionized 
gas all following the bar orbits. However, on the western 
side (mostly between —10 and +40"), the Nal CI follows a 
distinctly different, shallower velocity gradient. In pseudoslit 
2, the Nal CI gradient is flatter still, and between —20 to 
+20" is followed by the CO. These shallower gradients imply 
that this neutral and molecular gas must be located nearer to 
us along the line-of-sight (i.e. at larger galactocentric radii) 
where the radial component of the orbital velocity vector is 
less. 

Between —65 and 0" in the major axis pseudoslit we 
find a consistent area of blueshifted (10-60 kms -1 ) Nal C2 
points, and between +10 and +45" we find a cluster of red- 
shifted Nal C2 points (30-70 kms -1 ). Example line profiles 
and fits for these regions are shown in Fig. [2}} and c. In 
the velocity map (Fig. |4| we see that these redshifted C2 
velocities are present throughout the western section of the 
disk (as mentioned above). Together, these points appear to 
trace a component of neutral material with very different 
kinematics to the rest of the gas. 



3.2.1 Minor axis position-velocity graph 

Fig. [8] shows position-velocity plots from two 7"-wide 
pseudo-slits oriented along the minor axis (left panel) and a 
parallel line 20" to the east (pseudoslit 4; right panel). The 
location of these pseudoslits are also indicated in Figs.[4]and 
[6] As for the major axis plots, we include our new Nal and 
CO measurements together with the Ha emission line veloc- 
ities from |Paper I| We also include the Ha, [Nn]A6583 and 



the CO stays at a very constant velocity and is considerably 
redshifted from the ionized gas. In the north of pseudoslit 4, 
CO C2 traces Ha, whereas CO CI stays at a very constant 
velocity (from ~15-35"), remaining on the blue side of C2. 
The position of these points in the two p-v diagrams suggest 
that they may be tracing the far-side of the northern wind 
outflow cone, although their very constant velocities are sug- 
gestive of more quiescent gas clouds in the outer disk. We 
discuss this below in Section [4.21 

From —10" southwards the Nal and CO velocities agree 
well with ionized gas velocities, suggesting that here both 
the Nal and CO trace the far-side of the wind cone. At the 
equivalent offsets in pseudoslit 4, the CO emission is com- 
posed of two distinct velocity components. The redshifted 
component (C2, becoming CI) accelerates at a ~constant 
rate from —5" out to >— 35" (>600 pc), consistent with a 
number of Ha points between and —15". Between and 
— 15", the bluer component (CI; here blueshifted from C2 
by >100 kms -1 ), traces the kinematics of Nal CI and Ha 
, then jumps to ~— 110 kms -1 from —20" southwards (and 
assigned to C2). These split line components are consistent 
with tracing the wind outflow, as was noted by [Walter et al.| 
(20021. 



4 DISCUSSION 

M82 is a complex system that is not easily interpreted. The 
kinematics of any one gas phase in the central starburst re- 
gion are complex down to the resolution limit of any observa- 
tion made to date. Understanding the relationship between 
the different gas phases is therefore challenging, particularly 
when the sensitivity, spatial resolution and/or observation 
method of the data are not matched. In this study we have 
attempted to compare the neutral Nal absorption and 12 CO 
emission kinematics to the ionized (Ha) emission line kine- 
matics at the same spatial resolution, employing the same 
line fitting method for each. In the following we discuss some 
of the salient points from our results in context with the lit- 
erature and the M82 system as a whole. 



4.1 Neutral gas in the wind 

We find no evidence of Na I absorption that is kinematically 
associated with warm (~10 4 K) ionized gas outflow from our 
DensePak data. This could be due to (1) there being no cool 
neutral gas in the outflow, (2) a lack of velocity resolution 
meaning that weak blueshifted outflow components are un- 
resolved or washed out in our spectra, or (3) that the high 
inclination of the system means that the absorption line is 
confused with multiple line-of-sight components. Regarding 
point (1), it would be rather surprising that the outflow does 
not entrain any cool neutral material, given that we know 
there is ionized and molecular gas in the flow. 

We already know from Nal radial velocity measure- 



ments along sight-lines towards star clusters (Konstantopou- 



los et al. 2009 1 that the Na I is decoupled from the stars, 



clusters, and ionized gas. Although the clusters measured 



by Konstantopoulos et al. ( 2009 I were not strictly along 



Sra]A9532 data from |McKeith et al.| (l995) 



From +10 northwards along the minor axis (left panel) 



the major axis, the overall Nal radial velocity distribution 
within the central ~2 kpc of the disk is consistent with a 
flat gradient, suggesting that the absorbing neutral material 
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is located further out in the disk (where the radial compo- 
nent of its orbital velocity is smaller) . The high wind speeds 
in M82 should rapidly fragment clouds, so perhaps it would 
not be a surprise for Nal absorption from entrained clouds 
to be patchy (this effect can be seen in simulations such as 
Cooper et al. 1120091). Having said that, strong signatures of 



entrained Nal absorbers are found in many other starburst 
systems (e.g. |Rupke et al.|2005||Chen et al.|2010 l. 

On points (2) and (3), the velocity resolution of our 
DensePak spectra (~45 kms' 1 ) certainly does prevent us 
from being able to see the finer details of the Na I absorp- 
tion profile, known from previous higher resolution observa- 



tions ( [Smith k. Gallagher 2001] |Schwartz fc MarthT||2004 |. 
In Fig. [9j we compare the Nal line profile of M82-F, a su- 
per star cluster to the west of the nucleus (see Fig JT|), from 
our DensePak data to that measured by Smith & Gallagher 
(2001[) using the Utrecht Echelle Spectrograph (UES) on 



the William Herschel Telescope (WHT) at a resolution of 
8 kms' 1 (R=45000). Whereas five sharp absorption com- 
ponents can be identified in the UES spectrum, at the lower 
resolution of our DensePak spectrum only two (broader) 
components can be resolved. 



Schwartz fc Martin ( 2004 1 observed a ~4"xl'.'5 sight- 



line towards the centre of M82 clump A (centred at coor- 
dinates 09 h 55 m 52.6 s , +69° 40' 47'.'0, J2000; priv. comm.) at 
11 kins' 1 (R=30000) resolution with Keck/HIRES. They 
also found an extremely complex absorption system best fit- 
ted with five separate components, the two strongest being 
blueshifted by —90 and —35 kms -1 . 

These higher resolution observations show that the neu- 
tral gas kinematics in M82 are extremely difficult to dis- 
entangle, and great care must be taken. The presence of 
these narrow blueshifted components (in both the UES and 
HIRES spectra) suggests that there is indeed neutral gas 
entrained in the outflow which cannot be seen at the reso- 
lution of our DensePak spectra. Furthermore, the presence 
of so many components over a large range in velocities be- 
lies the fact that, due to the high inclination of the galaxy, 
these absorption profiles are likely to be quite confused by 
material along the line-of-sight, particularly in these central 
regions of the disk. 



4.2 Molecular gas in the wind 



By reanalysing the CO data from Walter et al. (20021 and 



performing a detailed multi-component line fit, we find con- 
clusively that molecular CO is entrained in the wind in these 
inner regions. In the southern wind plume (Fig. [8]), we iden- 
tify two components of CO emission, equivalent to the line 
splitting seen in |Walter e t al. (20021. The strongest compo- 
nent traces the redshifted (far-side) wall of the outflow cone, 
however a weaker component (C2) also traces the blueshifted 
side in some places (as seen in the pseudoslit 4 p-v plot). 
These results represent clear evidence that the inner wind 
contains entrained molecular gas that is coupled kinemati- 
cally to the ionized gas phase. 

A region of blueshifted (-100 kms' 1 ) CO emission can 
be seen in the south of the C2 velocity map (Fig. [6| at 
coordinates around (—30,-10), or in the pseudoslit 4 p-v 
plot (Fig. [8] right panel) between offsets of —20 and —35. 
By comparing this to the minor axis p-v diagram shown by 
Walter et ai] ( |2002[ ), we see that this component is equiva- 



lent to the faint blueshifted "prong" of the p-v fork on the 
southern side (their fig. 4). We therefore associate this com- 
ponent with the wind. Since this "prong" is only seen at low 
surface brightness, we do not identify it as clearly here since 



we masked low S/N data in the CO cube (see Section 2.2 l. It 



is interesting to note that this component is not seen directly 
along the minor axis, but only to the east. 

On the opposite (north) side of the disk, we also find 
redshifted (30-100 kms' 1 ) CO emission in both CI and C2 
(seen in both minor axis pseudoslits). These velocities are 
partly consistent with the ionized gas measurements (par- 



ticularly those from Paper 1 1, although the measurements 



of McKeith et al. ( 1995 I and Shopbell & Bland-Hawthorn 



( 1998 1 that extend further out indicate a much faster accel- 



eration in the cone back-wall than is seen in CO (see fig. 4 
of Walter et al.|2002 |. The fact that these CO components 
stay at very constant velocity with increasing disk height 
implies that they may trace more static clouds in the inner 
halo that could be pre-existing or tidal in origin. 



Wills et al. ( 2002 1 found evidence for four expanding 



shells in the inner disk of M82 from their high-resolution 
VLA A-configuration Hi data. The CO velocity splitting 
seen on the eastern side of the disk (see the major axis p-v 
plots in Fig. [7| is consistent with that measured by |Wills| 
et al. ( 2002[ ) for their shells 1 and 2. On the western side, 
however, CO counterparts to the Hi shells 3 and 4 can- 
not be identified, although CO line splitting consistent with 
the known molecular superbubble ~10" west of the nucleus 
( Weifi et al.|1999[ ) is clear. 



4.3 A counter-rotating Nal component 

Throughout the central regions (IFU positions 1,2,4), we see 
what appears to be a counter-rotating neutral gas compo- 
nent in the form of Na I C2 (see profile+fit examples Fig. [2)3 
and c). These Nal C2 velocities can clearly be seen in the 
major axis p-v diagram (Fig. [7]), and do not match anything 
seen in the ionized or molecular gas. Furthermore, the extent 
of these apparently counter-rotating Nal C2 velocities do 
not correspond to anything seen in optical imaging, includ- 
ing the dust lanes seen in the HST I-band image (Fig.0, 
or the CO map. 

M82-F is located at a projected distance of ~30" to 
the west of the nucleus and slightly above mid plane, just 
within the region where we find these unusual Na I velocities. 
As mentioned above, |Smith fc G allagher (2001) observed 
M82-F at high spectral resolution (R=45000; 8 kms" 1 ) and 
found a number of discrete narrow components in the Nal 
profile. Two of them are redshifted at 15 and 45 kms' 1 , and 
these correspond well with the C2 component identified in 
our DensePak spectra, which in this location is redshifted 
by 20-40 kms' 1 (as sh own in Fig. Ejl. 

(120001) presented V:A x I". 2 spatial resolu- 



Wills et al. 



tion VLA observations of M82 with a velocity resolution of 
~10 kms' 1 . For the few SN remnants (SNRs) that were 
clearly separated from other nearby sources, they were able 
to detect multiple Hi absorption components and measure 
their velocities. SNR 39.10+57.3, located ~ 5" west of the 
nucleus, 5" below the mid-plane, was found to have four ve- 
locity components at +60, +35, —10, and —120 kms' 1 (rel- 
ative to systemic). These redshifted components are consis- 
tent with our Nal C2 velocities. The Hi profiles against the 
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another three SNR discussed in their paper located ~25" to 
the west of nucleus (41.95+57.5), and ~20" and ~25" east 
of the nucleus (45.17+61.2 and 6.70+67.1), however, did not 
exhibit any absorption components outside the allowed val- 
ues of disk rotation. 

What could this counter-rotating component represent? 



eith et al.|1995 Paper 1 1. In Nai, however, we find no clear 
evidence of an equivalent outflow component. 



Smith & Gallagher (2001 \ comment that "this gas could be 
associated with infalling tidal debris or gas outflow." Al- 
though the outflow interpretation is difficult to reconcile 
with the other known wind tracers, the possibility that this 
component represents gas inflow is very interesting. Some 
simulations of galactic-shock accretion of gas in an oval po- 
tential, such as a nuclear bar, have produced highly eccen- 
tric streaming motions toward the nucleus, some portion of 
which are counter-rotating ( |Wada et al.|1998| ). 

The H I gas filaments in the halo of M82 are found i n the 
velocity range of —70 to 70 kms -1 relative to systemic ( Yun, 
Ho fc Lo]|1993 1 . Some of the Nai features could, therefore, 
originate in extraplanar gas. In particular the anomalous 
Nai velocities on the western side of M82 could represent 
an extension of the north-western Hi filament across M82 
( Yun et al.|1993 1. As the south-western Hi filament appears 
to join M82 at velocities close to the rotation speed, it is 
difficult to unambiguously distinguish between absorption 
from this filament and the outer disk of M82. The possibil- 
ity remains, however, that the south-western H I filament is 
feeding gas into M82. 

Nevertheless, the edge on inclination and the superpo- 
sition of disk and bar orbits with complex projection effects 
make it impossible to be sure from these data whether this 
counter-rotating component is indeed a sign of gas inflow 
or not. Higher sensitivity Hi imaging at good spatial and 
velocity resolution of the central few kpc of the disk would 
be needed to investigate this further. 



4.4 Consequences for interpreting Nai outflows 
vs. ionized gas 

Gas outflows have been identified in many star forming 
galaxies from low to high redshift (jHeckman et al.|1990||Pet-| 



tini et al. 2002 Tremonti et al. 2007 Weiner et al. 2009) 



Because of observational constraints it is common to mea- 
sure a single outflow speed for a galaxy from either an in- 
tegrated spectrum or a single long-slit pointing. Typically 
these are derived from measurements of interstellar absorp- 
tion lines (such as Nai, e.g. Rupke et al. 2005. Chen et al.| 
20101, but it is also common to use ionized emission lines 
(e.g. Forster Schreiber et al.|2009 i. Very few studies discuss 
the implications of using neutral vs. ionized and absorption 
vs. emission lines on determining the presence of galaxy out- 
flows and their velocities. 

This study highlights the inherent complexities of the 
relationship between the ionized and neutral phases, and 
therefore the assumptions of similarity between these phases 
and techniques. M82 is one of the nearest and best studied 
examples of a starburst and superwind system, and one that 
is often used as a template for more distant objects that can- 
not be resolved as well. However, spatially resolved Nai and 
Ha kinematics show that they cannot be used as equivalent 
tracers of the same gas. In Ha, a double-peaked line profile 
characteristic of the expanding superwind is clearly detected 
10-20" (few 100 pc) above and below the midplane (|McK 



In Fig. 10 we show the summed Nai and Ha- [Nil] pro- 
files over DensePak IFU positions 1, 2 and 4 (i.e. the cen- 
tral ~1 x 1 kpc of the disk), together with our best fitting 
multi- Gaussian fit and the fit residuals. For Nai, only one 
absorption component can be identified, with a velocity cen- 
troid of +23 kms -1 (relative to systemic) and FWHM = 
174 kms" 1 (mostly reflecting broadening due to rotation). 
Since no blueshifted component is seen, from this spectrum 
alone one would conclude that M82 does not host an out- 
flow. In contrast, three line components can be identified in 
the Ha and [Nil] lines, including two separated narrow com- 
ponents characteristic of outflowing material, and a broad 
(FWHM = 320 kms" 1 ) underlying component. 

We have shown that at higher resolution the absorption 
profile can break up into many narrow components (Sec- 
tion 4.1 1, and it may be that the summed profile would do 



the same. However, not all studies looking for absorption 
line outflows are done at high resolution (for example the 
DEEP2 spectra used by [Weiner et al. 2009 and Korne i et al.| 
|2012| were taken at R=5000, which is similar to this study). 
Even then, the lines-of-sight that can be probed by absorp- 
tion line studies rely on a bright background source, and this 
may not correspond with a region in which an outflow can 
be detected. For example, the high-resolution pencil-beam 



studies of Smith & Gallagher ( 2001 1 and Schwartz & Martin 
(2004} did not cover regions in M82 where the Ha outflow is 



most obvious (|Shopbell fc Bland-Hawth orn 1998) |McKeith 
et al.|1995||Westmoquette et al.||2009a i 



In conclusion, the high spectral resolution (R> 10000) 
observations are clearly preferable in order to search for sig- 
natures of a wind, so that the Na I line profile is sampled at 
a sufficient level. We have shown that high velocity compo- 
nents are hard to find in highly inclined systems, and may 
easily be missed in observations of similar systems at a sim- 
ilar spectral resolution to ours. 



5 SUMMARY 

We present spatially resolved (3" resolution) observations of 
the interstellar Nai D AA5890,5896 doublet absorption ob- 
tained with the DensePak IFU on the WIYN telescope, and 
compare them to the 12 CO (l->0) (2.6 mm, 115 GHz) emis- 
sion obtained with the OVRO millimetre array, across the 
inner ~2 x 1 kpc of the disk of M82. By measuring the Nai, 



CO, and ionized gas (Ha) kinematics (from Paper 1 1 at the 
same spatial resolution, and employing the same line fitting 
method, we have been able to make meaningful comparisons 
between the three gas phases. Our main conclusions are the 
following. 

In general, the Nai kinematics do not match those of 
Ha or CO, although in the eastern part of the inner bar 
region the three tracers do match one another closely, indi- 
cating that we are seeing molecular, neutral and ionized gas 
all following the bar orbits. The Na I rotation curve is flatter 
both along the mid-plane and at disk heights of a few hun- 
dred pc. The likely explanation of this is that the Nai line, 
being in absorption, probes gas that is located in front of the 
continuum source (the starburst), and is therefore biased to- 
wards material further out along the line-of-sight and thus 



© 2012 RAS, MNRAS 000,[T|fT6l 



Kinematics of the multi-phase interstellar medium in M82 7 



with a lower radial component of its orbital velocity. Despite 
this, however, the rotation seen in Nal CI is consistent with 
same offset axis as that found from the radio and optical 
recombination lines ( Rodriguez- Rico et al.|2004}|Paper I) of 
PA=12°. 

We detect a component of the Nal line with veloci- 
ties that are forbidden by the known orbital rotation of the 
galactic disk. Hints of this component were seen in high spec- 



tral resolution Nal observations of the cluster M82-F (Smith 
& Gallagher 2001 1 and the Hi absorption profile towards a 



SNR to the south-west of the nucleus ( |Wills et al.|2000| >. We 
interpret this as material that is counter-rotating or perhaps 
inflowing towards the starburst. These velocities have not 
been found in any observation at any other gas temperature 
phase, suggesting that this material must be cold. To inves- 
tigate this further, higher sensitivity H I observations of the 
inner few kpc of the disk are needed. 

We find a clear indication of entrained CO gas in the 
outflow with kinematics consistent with that of the Hq- 
emitting wind, thereby confirming the result of |Walter et al.| 



(2002|. It is possible that this CO line splitting follows the 



expanding shells detected in H I absorption measurements by 



Wills et al. ( 2002 1, and certainly traces the known molecular 
superbubble identified by ( |Weifi et al.| |T999). The existence 
of such structures might be expected to result from the mas- 
sive and highly clumped young complexes in the M82 star- 
burst zone. The breakout of hot gas from such high pressure 
bubbles would offer a likely source for the hot superwind. 

However, we find no evidence that Nal absorption is 
kinematically associated with the outflow. Since it would 
be rather surprising that the outflow does not entrain any 
neutral material (given that we know there is ionized and 
molecular gas in the flow), we conclude that a combination 
of (1) lack of velocity resolution and (2) confusion of the 
absorption profile due to the high inclination of the system 
is acting to prevent the detection. This view is supported by 
the observation of narrow, blueshifted Nal components at 
higher resolution (R>30000) along select sight-lines towards 



M82 ( |Smith fe Gall agher 200T] |Schwartz fc Martin|2004 1. 

Thus, in the search for neutral outflows from galaxies, 
high velocity components may easily be missed in observa- 
tions at similar spectral resolution to ours. This is particu- 
larly true for highly inclined systems. Often, the inclination 
of a galaxy is not known (especially for more distant sys- 
tems) so a non-detection of high velocity Nal absorption 
can simply be an inclination effect. Higher spectral reso- 
lution (R> 10000) is always preferable when studying ISM 
absorption lines, and particularly so in these situations. 
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Figure 1. HST / ACS F814W (I-band) image with the DensePa k IFU footprints o verlaid, together with arrows labelling certain well- 
known regions within the central starburst. The 2.2 um nu cleus (Lester et al,[ 



by a white cross. Contours show the integrated 12 CO flux ( [Walter et al. |2002} 



1990 



a = 09 h 51 m 43?4, 6 = 69°55'00", B1950) is marked 
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Figure 2. Examples of fits to the Nal lines, including the He I 5876A and telluric emission lines as labelled in panel (a). Zero velocity 
corresponds to the rest wavelength of Nal A5890. The y-ax.es are in arbitrary but relative flux units. The location from which these 
spectra originate in the IFU fields are indicated with the corresponding letters in the left panel of Fig. [4] 
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Figure 3. Map of the equivalent width of Nal absorption for the two line components. The (0,0) origin is centred on the 2.2 nm nucleus. 
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Figure 4. Nal line centroid velocity maps for components CI and C2 relative to v sys (=4-200 kms' 1 ). The (0,0) origin is centred on 
the nucleus. Letters in the left panel show the location of spectra shown in Fig. [2] The galaxy major and minor axes are marked with 
solid black lines. Pseudoslits along these axes, together with two additional pseudoslits (labelled slit 2 and 4) offset by 15" to the south 
of the major axis and 20" to the east of the minor axis (drawn in red), were used to make the position- velocity diagrams of Figs. [7] and 
[s] The pseudoslit widths are indicated by the dotted lines. The dashed line in the left panel represents the PA of the gaseous (Ha, [S III], 
P9) rotation axis (offset by ~12°; determined from the emission line maps in |Paper7| | . 
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Figure 5. Examples of fits to the CO line from selected locations as indicated with the corresponding letters in the left panel of Fig. [6] 
The fluxes are in units of Jy beam -1 . 



© 2012 RAS, MNRAS 000,[T|fT6l 



Kinematics of the multi-phase interstellar medium in 



M82 13 



a 




-60 -40 -20 20 40 -60 -40 -20 20 40 



arcsec arcsec 

Figure 6. 12 CO (1— >0) centroid velocity maps for components CI and C2, relative to v ays = 200 kms -1 . The overlaid contours represent 
the integrated CO flux, and clearly show the two molecular lobes. Letters in the left panel show the location from which the spectra 
shown in Fig. |5| were extracted. The same pseudoslits as shown in Fig.|4](used for the p-v diagrams of Figs. |7]andp} are shown. 
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Figure 7. Major axis position-velocity diagrams for Nal, CO and Ha l |Paper jj l, extracted from the corresponding 7" wide pseudoslits 
shown in Figs |4| and |6| Velocities are shown with respect to v ByB . The optical and near-IR ionized gas velocities measured along the major 
axis by |Cast lcs ct al. (1991) and McK eith et al.| l |1993| are plotted with + symbols. 
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Figure 8. Minor axis position-velocity diagrams for Nal, CO and Ha | |Paper Ip , extracted from the corresponding 7" wide pseudoslits 
shown in Fig s, |4| and |o] Velo cities are shown with respect to v aya . The optical and near-IR ionized gas velocities measured along the 
major axis by |McKeith et al.||l995| | are plotted with + symbols. The letter 'a' indicates a cluster of points that are discussed in the text. 
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Figure 9. A comparison of the Nal profile of M82-F, a super star cluster to the west of the nuclear regions (see Fig. ITJ, from our 
DensePak data (resolution 45 kms -1 ) and that measured by |Smith fc G allagher (2001) from their WHT/UES spectrum (resolution 
8 kms" 1 ). The velocity scale is relative to the blue component of the Nal doublet, and is shown relative to the systemic velocity of M82. 
The best fit profile and the individual components identified in the DensePak spectrum are shown in red. The vertical lines indicate the 
wavelengths corresponding to the systemic velocity for each of the Nal lines. Absorption components identified in the UES spectrum by 
|Smith &: G allagher (2001) are marked, together with absorption from the Milky Way (MW). In the DensePak spectrum, this Galactic 
absorption is filled by the telluric emission. 
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Figure 10. The summed Nal and Ha- [Nil] 

spectra from DensePak IFU positions 1, 2 and 4 (i.e. the central 
~lx 1 kpc of the disk), shown together with the best fit to the 
emission and absorption lines, the individual components to the 
fit, and the residuals. For Nal only one absorption component 
can be identified with a velocity ccntroid of +23 kms -1 (rela- 
tive to systemic) and FWHM = 174 kms -1 . In contrast, three 
components can be identified in the Ha and [N n] lines. 
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